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Abstract 

The Newman-Penrose formalism is used to deal with the massless scalar, neu- 
trino, electromagnetic, gravitino and gravitational quasinormal modes (QNMs) in 
Schwarzschild black holes in a united form. The quasinormal mode frequencies eval- 
uated by using the 3rd-order WKB potential approximation show that the boson 
perturbations and the fermion perturbations behave in a contrary way for the vari- 
ation of the oscillation frequencies with spin, while this is no longer true for the 
damping's, which variate with s in a same way both for boson and fermion pertur- 
bations. PACS: number(s): 04.70.Dy, 04.70.Bw, 97.60.Lf 

Keywords: Schwarschild black hole, QNMs, arbitrary spin fields, WKB approxi- 
mation. 



Ever since Chandrasekharfl] and Vishveshwara[2] discovered the quasi- 
normal modes of black holes, much effort has been devoted to investigating 
the QNMs of various black hole cases[3,4,5,6,7]. By obtaining quasinormal 
mode(QN) frequencies, we can not only test the stability of the spacetime 
against small perturbations, but also probe the parameters of black hole, such 
as its mass, charge, and angular momentum, and thus, help uniquely identify 
a black hole. 

QNMs are described as the pure tones of black hole. They are defined as so- 
lutions of the perturbation equations belonging to certain complex character- 
istic frequencies which satisfy the boundary conditions appropriate for purely 
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ingoing waves at the event horizon and purely outgoing waves at infinity [8]. 
QNMs are excited by the external perturbations(may be induced, for example, 
by the falling matter). They appear as damped oscillations described by the 
complex characteristic frequencies which are entirely fixed by the properties 
of background geometery, and independent of the initial perturbation. These 
frequencies can be detected by observing the gravitational wave signal [9] : this 
makes QNMs be of particular relevance in gravitational wave astronomy. 

QNMs were firstly used to study the stability of a black hole. Detweiler 
and Leaver found the relations between the parameters of the black hole and 
QNMs. Latest studies show that QNMs play an important role in the quantum 
theory of gravity. This is related to the quantization of black hole area[10]. 
For example, there exit some possible relations between the classical vibra- 
tions of black holes and various quantum aspects, such as the relation between 
the real part of the quasinormal mode frequencies and the Barbero-Immirzi 
parameter, a factor introduced by hand in order that loop quantum grav- 
ity reproduces correctly entropy of the black hole[ll,12,13]. All these works 
deal with asymptotically flat spacetimes. The recently proposed AdS/CFT 
correspondence makes the QNMs more appealing, due to its argument that 
string theory in anti-de Sitter(AdS) space is equivalent to conformal field the- 
ory(CFT)in one less dimension [14]. In addition to this context, many studies 
also have been done on QNMs of various spin- fields [3, 7, 8] . Chandrasekhar[8] 
investigated the QNMs of fields with spin s = 1/2,1,2 in Kerr black hole, 
and s = 2 in Schwarzschild and Reissner-Nordstrom black hole . Cardoso and 
Lemos[3] studied the QNMs of Schwarzschild- anti-de Sitter for fields with spin 
s = 1,2. Cho[7] calculated the massive Dirac quasinormal mode frequencies of 
Schwarzschild black hole. However, problem on how to deal with QNMs due 
to arbitrary spin field (s = 0, 1/2, 1,3/2,2) in a united form has never been 
discussed in these previous works . The main purpose of this article is to give 
a possible way to deal with QNMs in Schwarzschild black hole of arbitrary 
spin field in a united form. 

We start with the line element in standard coordinates for the Schwarzschild 
space-time 

ds 2 = -e 2U dt 2 + e~ 2U dr 2 + r 2 (d6 2 + stn 2 6d V 2 ) , (1) 

with 

e 2U = l , (2) 

r 

where M is the mass of the black hole. 

The Teukolsky's master equations[15,16] for massless arbitrary spin fields (s = 
0, 1/2, 1, 3/2, 2) in Newman-Penrose formalism can be written as[17] 

{[D - (2s - l)e + e* - 2sp - p*](A - 2s7 + //) 
-[5 - (2s - 1)P - a* - 2sr + n*](5 - 2sa + tt) 

_( s _l)( 2s -l)^ 2 }$ s = 0, (3) 
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and 



{[A + (2s - 1) 7 - 7* + 2s/i + fi*](D + 2se - p) 
-[8+ (2s - l)a + 13* + 2stt - r*] (5 + 2sf3 - r) 

-(s-l)(2s-l)tt 2 }$_ 8 = 0. (4) 

Assume that the wave-functions in Eqs. (3) and (4) have a t- and a in- 
dependence specified in the form e l ^ t + mL P) ) i. e ., 

$ s = R +s (r)A +s (0y {uJt+mip) , *-s = R- s (r)A_ s (9)e-^ t+m *\ (5) 
and define 

P +s = A S R +S , P_ s = r 2s R. S} (6) 

(where R± s and A± s are, respectively, functions of r and 9 only, and A = 
r 2 — 2M) one can decouple equations (3) and (4) as two pairs of equations, 



[A^s^ - 2 (2s - l)iuor]P +s = XP +S , (7) 
^fi_ s ^ s A +s = —XA +S , (8) 



and 



[A@l_ 8 @o + 2(2s - l)iur]P- s = AP_ S , (9) 
^ S ^A_ S ^-XA_ S , (10) 

where A is a separation constant. The reason we have not distinguished the 
separation constants which derived from Eqs. (3) and (4) is that A is a pa- 
rameter that is to be determined by the fact that A +s should be regular at 
9 = and 9 = ir, and thus the operator acting on A_ s on the left-hand side 
of Eq.(10) is the same as the one on A +s in Eq.(8) if we replace 9 by n — 9. 
In Schwarzschild black hole, the separation constant can be determined analytically [16, 18] 
for boson 



A = (1 + \s\)(l- \s\ + 1), 1 = |s|,| s | + l,... (11) 

for fermion 



A = 0' + |s|)0'-|s| + l), j = |s|, ]s| + 1, • • • , and j = l±\s\ (12) 

where / and j represent angular quantum number and total quantum number, 
respectively. Since P +s and P_ s satisfy complex-conjugate equations (7) and 
(9), it will suffice to consider the equation (7) only. 

By introducing a tortoise coordinate transformation dr\ = ^dr, one can 
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rewrite the operators @ and 3)\ as 



^2 ^,2 

-A+, and f? f = -A_. 



± iu;. 



(13) 



where we have defined A± = , 
With the definition Y = r 1 ^ 2s P +s , equation (7) can be written as 



„2s+3 



A 

where A 2 
the form 

where 
and 



A 2 Y + 



d , r 

dr* A 



Is 



A.y + 



d 



1 d 



,2s-l 



dr \ A s 1 dr 



- Ar 



2s-l 



d 2 



Y = 0, 
(14) 



d(r ^ + or. On further simplification, Eq.(14) can be brought to 

A 2 r + pa_f - gr = o, (15) 



„ d r M 
dr* A s 



<2 = ~2 



(16) 
(17) 



2A A' 

A -(2 S -1)( S -!)(—-- 

Considering the purpose of this article, we should seek to transform Eq.(15) 
to a one-dimensional wave-equation of the form 



A 2 Z = VZ, 



(18) 



where V represents potential. 

The transformation theory introduced in Ref. [8] is applicable to solve this 
problem. We assume that Y is related to Z in the manner 



Y = £A+A+Z + WA + Z, 



(19) 



where £ and W are certain functions of r* to be determined. One can then 
deduce the following equations [8] 



*r [ X ,X =^(QT-2^) + ft (21) 



where if is a constant, and Xi T are certain functions of r*. 

The following work is to look for solutions of equations (20)-(23). These 
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equations provide four equations for five functions £, /3, \i T, and V. As 
a result, there is considerable difficulty in seeking useful solutions of these 
equations. An obvious fact is that x an d V are independent of u (i.e., they 
do not contain any term linear in iuj ). Under these considerations, we can, 
without loss of generality, suppose that T, f3, K are of the forms [8] 

T = T 1 (n) + 2iujf(s), = p 1 (n) + 2iu J p 2 , K = K 1 + 2iuK 2 , (24) 

where (3 2 , K ii K 2 are constants and f(s) is function of s. In this article, we take 
the choice 

/( s ) = i s (2 S -l)(6 S 2 -23s + 23) far s = 0,^,1, | 2. (25) 

Making use of the fact that, for a equation contains real and imaginary parts, 
the real parts and imaginary parts in two sides of the equation are equal, 
respectively, we can separate each of the Eqs.(21)-(22) into two equations by 
substituting Eq.(24), i.e., 

X = fQ + ^fo, (26) 

A ( r As \ r As 

= (27) 

and 



/3i/ + /^7i = «2, (28) 

* ~ XjF, + = ^ (29) 

where k = K\ + 4u 2 f(32- 

Substituting Eqs.(26) and (28) into Eq.(27), we obtain 

Tl = (30) 

Here we have defined F = L ^-, and ',r*' denotes the differential with respect 
to TV Eq.(29) can then be written as 

^UF + h) -f fF _ (32 + {fF _p 2? (31) 

It's obvious that Eq.(31) is a condition on F if solutions of the chosen form 
are to exist, and hence the work to seek available [3 2 , k, and k 2 whose values 
satisfy Eq.(31) is a key step to obtain the function of potential V. Since n 2 
occurs as n 2 in Eq.(31), two choices which associated with +k 2 and —k 2 are 
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possible to satisfy the equation. Further study finds, an available choice is the 
following 



/3 2 = --(2s-3)(s-2)(4 S -l)A, 



k = -(2s - l)(s - l)(2s - 3)(5s - 8)A(A + s), 



K 2 = s(s - 1) 



-(2 - s)\Js - - + A + (2s - l)(2s - 3)M 



(32) 



The solution for V can then be obtained by substituting the expressions of (3 2 , 
k, and k 2 into Eq.(23),i.e., 



„4s 



(/^-/9 2 ) J 



(fF-p 2 r 



(33) 



where we have distinguished the transformations associated with +k 2 and — k 2 
by superscripts (±). 

It is obvious that y(+) equals to y* - - 1 for s = 0, 1 from the equation (33). 
For the case of s — |, |, 2, one can obtain a simpler form of yW by defining 
a new function F 

= ±k 2 ^ + k 2 2 F 2 + kF. (34) 



The expression of this new function is 



„4|s-l| 



A + 2(2s- l)(s- 1) 



M 



(35) 



References [8, 19] has shown that potentials and y( _ ) related in the way 
Eq.(34) shows are equivalent and hence possess the same spectra of quasinor- 
mal mode frequencies. We shall therefore concentrate on only in evalu- 
ating the quasinormal mode frequencies. 

So we can simplify the radial equation (7) to a one-dimensional wave- 
equation of the form 

+ u 2 Z = VZ, 36 

dri 

with smooth real potentials, independent of cu, i.e., 



V = —F 

r -4s 



( f F-(3 2 )fF^-pF^ 2 + K 2 F^ 

(fF - p 2 y 



(37) 



where F is a known function of r*. Note that we have written y(+) as V 
because we shall not work with which will give the same quasinormal 

mode frequencies. 

Figure. 1 demonstrates the variation of the effective potential V with spin 
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s. From this we can see that peak values of the effective potential V decrease 
with s for boson perturbations, while they increase with s for fermion per- 
turbations. This phenomena is closely related to the value of the separation 
constant A. 

Analytic expressions for the quasinormal mode frequencies are usually 
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2 2.5 3 3.5 4 4.5 5 Fig. 1: Variation of the effective 
r potential V with spin s. 

very difficult to obtain. One hence should appeal to approximation schemes 
to evaluate these frequencies. Many methods are available for our purpose. 
One often used is WKB approximation, a numerical method first proposed 
by Mashhoon[20], devised by Schutz and Will [21], and was subsequently ex- 
tended to higher orders in[22,23]. An obvious character of this scheme is that 
it is very accurate for low-lying modes (n < Z) [24] . Under these considera- 
tions, we may use WKB approximation to evaluate the quasinormal mode 
frequencies within low-lying modes. The values for n < I are listed in Ta- 
bles 1 and 2, which also include a comparison with the results of Cho[7], 
and those of Chandrasekhar and Detweiler(CD)[l]. Some gravitational modes 
(I — 3, n — 2, and I — A,n — 3) of CD's results have no counterpart in our 
results. This anomalous values were presumably results of numerical instabil- 
ities in their computational method that appeared when Im(o;)~ Re(u;) (see 
Ref. [1] for details). From Table 1, we see that our results for boson perturba- 
tions are almost the same as Iyer's[24]. Our values for Dirac modes agree well 
also with the results of Cho from Table 2. Notice that during our evaluating 
procedures, we have let the mass M of the black hole as a unit of mass so as 
to simplify the calculation. For negative n, the values are related to these with 
positive n by reflection of the imaginary axis. Figure 2 and figure 3 depict the 
variation of the QN frequencies with spin s. Some significant results can be 
found from these figures, (i) For boson perturbations, the real parts of the 
QN frequencies decrease with the spin value for the same I and n, while they 
behave in a contrary way for the fermion perturbations, (ii) The imaginary 
parts of the frequencies decrease slowly with spin value, regardless of boson 
or fermion perturbations, as showed in Fig.3. This means that the variation 
of the oscillation frequencies with spin for boson perturbations is totally dif- 
ferent from fermion perturbations, while the variation of the damping's with 
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spin has a same law for both boson and fermion perturbations. 

The method we used to deal with QNMs of arbitrary spin fields in this ar- 
ticle is also available for discussion on other black hole cases, such as Reissner- 
Nordstrom and Kerr black hole, but a united form to deal with arbitrary spin 
perturbations in these black holes is still an unsolved problem. 
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Table 1 

Quasinormal mode frequencies for boson perturbations. u>i(i = 0,1,2) represent 
scalar, electromagnetic and gravitational perturbations, respectively, locd represent 
Chandrasekhar's results for gravitational perturbations. 
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n 




^0 



















1046+0 


11521 








1 








2911+0 


0980i 


0.2459+0.09311 






2 








4832+0 


0968i 


0.4571+0.09511 


0.3730+0.08911 


0.3737+0.0890i 




1 





4632+0 


2958i 


0.4358+0.29101 


0.3452+0.27461 


0.3484+0.2747i 


3 








6752+0 


0965i 


0.6567+0.0956i 


0.5993+0.0927i 


0.5994+0.0927i 




1 





6604+0 


2923i 


0.6415+0.2898i 


0.5824+0.2814i 


0.5820+0. 2812i 




2 





6348+0 


4941i 


0.6151+0.4901i 


0.5532+0.4767i 


0.4263+0.3727i 


4 








8673+0 


0964i 


0.8530+0.09591 


0.8091+0.09421 


0.8092+0.09411 




1 





8557+0 


2909i 


0.8411+0.28931 


0.7965+0.28441 


0.7965+0.28441 




2 





8345+0 


4895i 


0.8196+0.48701 


0.7736+0.4790i 


0.5061+0.4232i 




3 





8064+0 


69261 


0.7909+0.68921 


0.7433+0.67831 





Table 2 

Quasinormal mode frequencies for fermion perturbations. = 1/2,3/2) represent 
Dirac, Rarita-Schwinger perturbations, respectively, ujcho represent Clio's results 
for Dirac perturbations. Notice that k = I + 1 for j = I + 1/2, according to Cho's 
definition. 
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n 




^1/2 


^Cho 


^3/2 
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2 








3786+0.09651 


0.379+0.0971 
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3 








5737+0.09631 


0.574+0.0961 


0.7346+0.0949i 






1 





5562+0.29301 


0.556+0.2931 


0.7206+0.28701 
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4 








7672+0.09631 


0.767+0.0961 


0.9343+0.0954i 






1 





7540+0.2910i 


0.754+0.29H 


0.9233+0.2876i 






2 





7304+0.4909i 


0.730+0.49H 


0.9031+0.4835i 


4 


5 








9602+0.0963i 


0.960+0.096i 


1.1315+0.0956i 






1 





9496+0.2902i 


0.950+0.2901 


1.1224+0.2879i 






2 





9300+0.48761 


0.930+0.4881 


1.1053+0.48281 






3 





9036+0.68921 


0.904+0.6891 


1.0817+0.68121 
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Fig. 2. Variation of the QN frequencies with spin s, the left for boson perturbations, 
the right for fermion perturbations. 
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Fig. 3: Variation of the imaginary 
parts of the QN frequencies with 
spin for 1 = 4. 
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